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A state of a black hole in loop quantum gravity is given by a distribution of spins on punctures 
on the horizon. The distribution is of the Boltzmann type, with the area playing the role of the 
energy. In investigations where the total area was kept approximately constant, there was a kind of 
thermal equilibrium between the spins which have the same analogue temperature and the entropy 
was proportional to the area. If the area is precisely fixed, however, multiple constraints appear, 
different spins have different analogue temperatures and the entropy is not strictly linear in the 
area, but is bounded by a linear rise. 



I. INTRODUCTION 

Black holes are generally regarded as pure manifestations of gravity. Apart from the obvious geometric properties of 
black holes, they have been known to show intriguing thermodynamic features even though at the classical level they 
do not allow anything to escape and cannot have non- vanishing temperatures. The area of the horizon of a black hole 
was shown to behave like an entropy in [1| . Subsequently the quantum field theory of a particle in the field of a classical 
black hole led to the assignment of a temperature to a black hole and the area was quantitatively interpreted as a 
measure of the thermodynamical entropy. There have been many attempts to understand this entropy in a quantum 
theory of gravity. Recently, the framework known as loop quantum gravity has yielded a detailed prescription for 
counting of microscopic quantum states corresponding to a black hole The quantum states are associated with 
cross sections of the horizon carrying some punctures. Spin quantum numbers j, m attached to the punctures label the 
quantum states. The entropy is obtained by directly counting the possibilities of labels that are consistent with a given 
area and has been seen to be proportional to the area in different approaches (1, 0, [E, Q . Numerical counting of the 
number of states however has exhibited a modulation of this linear behaviour [7| . This is apart from the logarithmic 
corrections that are also seen but are well understood 0, d, Q . The breakdown of strict linearity indicates that there 
is something beyond the calculations of 0, H, @] ■ Some attempts have been made to explain this departure from 
linearity [9|, il^] . But a fundamental issue is involved in connection with the area. Whereas the earlier predictions of a 
linear rise in the entropy were made by keeping the area only approximately fixed, the numerical studies fixed the area 
very precisely. This necessitates a fresh consideration of the area constraint. The result is an entropy that depends 
on the area in a more complicated way where it is bounded by a linear rise and saturates the bound at isolated points. 

In the analyses of [1, S @, where the area was only approximately fixed, the irrational nature of the quantities 
\/J(jr^-l) appearing in the eigenvalues of the area operator was ignored: essentially, j, which is a half-integer or 
integer, was taken for convenience as just a real variable. In the present calculation, where the area is precisely 
fixed, the irrational nature of each + 1) must be taken into account. This necessitates a classification of spins 

into isolated classes of compatible spins as explained below. Spins of each class are in a kind of thermal contact and 
equilibrium: the distribution of the total area amomg different spins is governed by analogues of temperatures, one for 
each class. The effective analogue temperature depends on the ratios of the contributions to the area of the different 
classes of spins. Its variation can explain the oscillations in the variation of the number of states with the area. 

We first do the counting treating the quantum numbers j as well as the m as labels characterizing quantum states, 
as in [5|. This is our preferred count, because the punctures on the horizon are assigned both quantum numbers j 
and m [51]. Thereafter we count the number of states labelled by only the m, which is a popular counting criterion 
relying only on the surface Hilbert space 0, @]- We end with various concluding remarks. 

II. COUNTING OF STATES LABELLED BY QUANTUM NUMBERS j, m 

It is convenient to use units such that Anji'j, = 1, where 7 is the so-called Barbero-Immirzi parameter involved in 
the quantization and £p the Planck length. Setting the classical area A of the horizon equal to an eigenvalue of the 
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area operator for a specific spin configuration of punctures on the horizon, we write 

A = 2^s,-,„V7(JTT), (1) 

where Sj^m is the number of punctures carrying spin quantum numbers j, m. Such a spin configuration is admissible 
if it obeys ([T]) together with the spin projection constraint 



msj^jn, (2) 



which is required for a quantum horizon with the topology of a 2-sphere; it arises from the consistency of a bulk Hilbert 
space and a surface Hilbert space in the theory [3] . The total number of quantum states for these configurations is 

ds,.^ = n , , ■ (3) 

To obtain the dominant permissible configuration that contributes the largest number of quantum states, one may 

maximize Indsj^ by varying sj^m subject to the constraints. In the area constraint and the spin projection were 

considered. Now we have to understand that the strict area constraint is a very severe one and may be decomposed 

into several constraints. This is because in varying the states, we can vary the Sj^m only integrally and a sum of 

these integers multiplied by the irrational factors y/ j{j + 1) has to be kept unchanged. Not all j-s can mix with one 

another in such a variation. The quantities \/ j{j + 1) for different j are in general relatively irrational, i.e., have 
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irrational ratios, except in special cases. For instance, the quantities for j = ^, i = 3 j = ^ and j = 48 are in the 



ratio 1:4:15:56; again, the quantities for j — 1, j — 8 and j — 49 are in the ratio 1:6:35, the quantities for j — ^ and 
j = 15 are in the ratio 1:8 and the quantities for j = 2 and j = 24 are in the ratio 1:10. Thus the set of values of 
j gets divided in to disjoi nt subsets such that the j in each subset are compatible with one another in the sense that 
the quantities y^j{j + 1) are in rational ratios within a subset but in irrational ratios in different subsets. When we 
vary the Sj.m, mixing can occur only within such a subset but not across subsets. So the contribution of each subset 
to the total area must remain fixed: 

An = 2 Sj,m\/j{i + 1) const, (4) 

where N denotes such a subset. 
Now one can write 

5 In ds^^ = In ^ Sj,™ - ^(i5sj,m In Sj,m) (5) 

for small 5sj^rn- The condition for the maximum can be expressed in terms of Lagrange multipliers Aat, ol: 



In ^ -2\N^]{j + 1) - am, ( for 3 e N) (6) 



whence 



^LUL^ = g-2A„VjO+i)-"™ (for jeN). (7) 



Consistency requires that Xjy and a be related to each other by 



^^g-2A„V7aTT)^g-am^^_ (8) 

N jeN m 

In order that ([7]) satisfies the spin projection constraint, we need J2m "le"""^ = for each j, which essentially implies 
a = 0. Therefore, the consistency condition becomes 



J2 e-^^-VJU^ = l. (9) 

N jeN,m 
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Here the sum over N goes over the subsets appropriate for the area in question. 

This is similar to the equation obtained earher using a single area constraint. However, while that equation could 
be solved for its single A, now there are several variables Xn in general. There are also more equations: 



2(E ^^■.") E V7(JTT)(2j + l)e-2^«V^, (10) 
jeN 

where the 2j + 1 comes from summing over m. This implies 



Aj, E.eN VW+^m + l)e-^^"V^ ^^^^ 

^ En E.^n V7(7 + T)(2j + l)e-2^« ' 

which, together with ([9]), determine the Ajv- 
Note that in Stirling's approximation. 

In d.sj „^ = E "^nAn + a E '^j.™'™' (12) 

N 

in which the last term vanishes, so that the leading contribution to the entropy is 

^ = E^^r^ (13) 

if normal units are used. 

It has to be understood that the a-re determined by A and the Ajv are in turn determined by these. It is not 
possible to solve the equations explicitly in general, but let us consider some simple cases. The simplest possible case 
involves a single subset N with i = 5, J = 3, j = j = 48, .... This corresponds to the area being an integral multiple 
of \/3 in our units. Numerical solution of ([9]) for this special case yields 

A^^'"^ = 0.521..., (14) 

which is greater than « 0.400... corresponding to just spin 1/2 [s*], but is less than the 0.86 corresponding to the 
use of a single area constraint @. 

Next we take two A''-s, one with j = 51 i = 3, j = 48, ... and the other with j = 1, j = 8, j = 49, but keep 

only the smallest values of j in each subset as an approximation. This means that the area must be the sum of an 
integral multiple of a/S and an integral multiple of 2\/2 in our units. Then 

Ai oc \/3 exp(— Ai\/3) 

A2 (X 3V2exp(-2A2V2). (15) 

These equations, together with 

2e-Ai\/3^3g-2A2\/2 ^ (16) 

determine Ai, A2 as functions of Ai/A2. The entropy, involving these A-s, can be written as S* = XA with an average 
A. We find that as A1/A2 varies, A increases from and falls to y= after reaching 

X\max = 0.704 at A1/A2 = 0.884, (17) 

corresponding to Ai = A2. 

If several N are involved, finding the corresponding A at is on the same lines, but more complicated. The average 
A again varies with the ratios of the A^-s and can be seen to be maximum when all Xn are equal. The value of 
this maximum depends on the subsets N involved, and reaches its peak value only when all spins participate; then it 
becomes equal to 0.86 by virtue of Thus A is bounded by the value 0.86 corresponding to a single A. 

In general, as the area A is varied, the relevant subsets N change and so do the Xn- For values of A which are not 
sums of eigenvalues, there is no state at all, but when A can be expressed as a sum of terms of the form \/ j{j + 1), 
some Ans are nonzero. The number of states involves also the corresponding Xn- The quantities An,Xn fluctuate 
with variations of A and so does A, which is not a constant but depends on the ratios of the 'components' A^r of A. 
That is why the variation of the number of states is not monotonic and a structure is seen. As A < 0.86, the plot 
is generally below the S — 0-86A line, but touches it at points where all spins participate and moreover all Xn are 
equal, which happens for special ratios of the A^v-s. A falls where the area is not an eigenvalue. Such oscillations with 
a linear bound have been seen numerically 0. A step- like structure emerges on bunching of the area variable into 
bins because of the discrete nature of area eigenvalues. 
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III. COUNTING OF STATES LABELLED BY QUANTUM NUMBERS m 

The ideas of the previous calculation can be easily extended to do a fine-grained counting of states labelled by only 
the m quantum numbers as envisaged in . Here one has to consider Sm = j ^^'^ maximize the combinatorial 
factor by varying the Sj,m- Maximization of the combinatorial factor occurs not in the interior but on the boundary 
of the configuration space @. The consistency condition for the Lagrange multipliers Ajv becomes 



^^(2 + ,5,i)e-2^«V^(^ = l, (18) 



where for each m one needs the minimum j possible. Furthermore, 

-2\Ny/jU + l) 

(19) 



An _ E.eN V7(7 + T)(2 + ^,i)e-^^"V^ 
^ En E.eN vW + T)(2 + S,,)e-^^-VJuTI)' 



which, together with the preceding equation, determine the Ajy. 

If one has an A such that only one Ajy is nonzero, corresponding to the subset with j = — 3, j — — 48, 
the corresponding Ajv can be easily determined from the consistency condition to be 

A^™^ = 0.453... (20) 

This is less than both the value 0.521 obtained above with j, m counting and the value with a single constraint but m 
counting, namely 0.790 Q. The case of two nonzero An with j — and j = 1.. is the same as in the counting with 
j,m labels and again produces a peak A value of 0.704. When there are several nonzero An the number of states 
is determined by the corresponding Xn which can in principle be obtained from the above equations. The An,Xn 
fluctuate with change of A and produce a structure as in the previous situation, with A < 0.79. 



IV. CONCLUSION 



The earlier approximate analytical calculations, valid for large areas, indicated that the entropy, increases linearly 
with the area of the horizon up to a logarithmic correction. In this Letter a more precise analytical calculation has 
been carried out, leading to an understanding of a numerically observed modulation of the linearity. The point is 
that when one talks about the area eigenvalues, these involve irrational numbers \/ j{j + 1). Real numbers in general 
cannot be expressed as a combination of numbers of this sort, but some can be. Such numbers can often be expressed 
in many different ways as combinations of this form. The building units \/ j{j + 1 ) fall naturally into disjoint classes. 
Members of a single class have rational ratios, so that different combinations of members in a class can be constructed 
with a constant sum. Members of different classes however have irrational ratios. If the total area is to be kept 
constant, the contribution from each of these classes has to be kept constant separately. Thus there is a number of 
constraints corresponding to the different classes. As the number of punctures is finite for finite A, the number of 
subsets or constraints also must be finite. In practice the effective number will be small because only low spins j 
can contribute significantly to the counting of states. In any case, each class or each constraint is associated with a 
Lagrange multiplier which can be determined if the contributions of the different classes to the total area are known. 
The number of states can be expressed as a function of all these variables: the area contributions from the different 
classes and the Lagrange multipliers. While they depend on and are determined by the total area, they fluctuate as the 
total area cha nges: indeed the contributions of the different classes to the total area change wildly as different classes 
of \/ Jii + 1) involved. On top of this there are values of the area which cannot be expressed as combinations 
of V jU + 1) s-'^d so yield no quantum state. These fluctuations result in a non-monotonic change of the number of 
states with the total area bounded by a linear growth. These oscillations go over into steps on consideration of bands 
of area. 

In the above discussion we have considered only the leading area dependence, ignoring the logarithmic corrections 
found in earlier papers. Clearly, the introduction of new constraints has a role to play in the logarithmic corrections. 
The usual term (— ^InA) arises because of the presence of two constraints, viz., the area and spin projection con- 
straints. Each additional constraint will contribute a — ^InA. However, the number of new constraints is finite for 
finite area, so the number of subsets N is only finite. It is only of academic interest to calculate N and the logarithmic 
corrections for a general case. 

The above discussion may cause some confusion about the correct value of the A or of the Barbero-Immirzi pa- 
rameter. It has to noted that we have calculated the degeneracy of a fixed area eigenvalue. In contrast, in earlier 
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work, one calculated the number of area eigenstates whose corresponding area eigenvalues fell within a narrow band 
[A — €,A + e], where e/A <C 1, around a fixed area A. When there is a band, the area constraint is implemented only 
approximately: one can always find a rational number close to + 1) within the band and the additional An 
constraints do not arise. Therefore the earlier calculations were appropriate in a statistical sense. They were realistic 
in the sense that the larger an area is, the larger are the errors in fixing it, hence a band is more appropriate. For 
a fixed eigenvalue the number of degenerate states must be smaller than the total number of states calculated for 
a band. This is related to the increase in the number of constraints in the present calculation. So the smaller and 
varying values of the parameters obtained in this paper are not unexpected. The earlier values are reached as the 
peak values of A as it varies with varying ratios of the A^s. 

Each A behaves like an inverse temperature - but it is dual to the area rather than to the energy, in the sense that 
A ~ The appearance of several Xn corresponds to a many-temperature system. This means that the system 
as a whole is not in what may be thought of as an analogue of thermal contact; it has different subsystems each 
characterized by its own analogue temperature. This occurs because area may be exchanged between spins in one 
class, but not between different classes. Of course, these analogue temperatures are not related in any way to the 
actual Hawking temperature. 
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